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Abstract. We study noise and noise energy of a high-Tc dc SQUID fabricated on a 
high-efl, substrate whose conduction properties are given by transmission line physics. 
We show that transmission line resonances greatly enhance the noise. Remarkably, 
resistance asymmetry enhances these resonances even more. However, as the transfer 
function (^^^) scales the same way, the noise energy is reduced by asymmetry greatly 
enhancing the flexibility and performance of the SQUID. 



1. Introduction 

SQUIDs are versatile high-resolution magnetometers that find wide application in 
fields ranging from brain research, low-field MRI, geological prospecting, precision 
electrometry, detection of elementary particles and quantum computing. They are so 
powerful because of their large sensitivity and low intrinsic noise, leading to excellent 
signal-to-noise ratios that in some cases approach the limit set by quantum mechanics. 
High Tc SQUIDs offer the possibility to be operated at temperatures of liquid nitrogen, 
lowering cost and enhancing flexibility. Their electrodynamics is more involved than 
that of their low-Tc counterparts and optimization of signal-to-noise figures requires 
renewed attention. 

In a previous paper [1], we have analytically studied transmission line resonances 
in high Tc dc SQUIDS. Such resonances are exhibited in the characteristics of SQUIDs 
fabricated on substrates with high dielectric constant like strontium titanate. In [1], 
we analytically derived the SQUID power balance equation for both symmetric and 
asymmetric SQUIDs and investigated SQUID current- voltage liV), voltage-fiux V^($) 
and voltage modulation AV characteristics. In this paper, we analytically study the 
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effect of transmission line inductance on tfie noise cfiaracteristics of a dc SQUID. We 
will closely follow the methods used in [2]. 

The paper is organized as follows. In section 2, we set notations by describing the 
relevant circuit equations. In section 3, we study the effect of asymmetry on white noise 
in transmission line dc SQUIDs. Wc conclude in section 4. The calculational details are 
provided in appendix A. The SQUID parameters are the same as those used in [1] . 



(2) 



2. Circuit equations 

The geometry of the SQUID washer that wc have used is shown in fig. (la) [1]. This is a 
geometry commonly used to manufacture SQUIDs [3, 4, 5] and is the geometry Enpuku 
et al used to numerically investigate the effects of large dielectric constant of strontium 
titanate (STO) on the characteristics of high Tc dc SQUIDs [6]. The slit of the SQUID 
washer makes up the SQUID inductance where /, s, w and d denote the slit length, slit 
width, electrode width and thickness of electrode respectively. For this geometry, the 
inductance per unit length of the slit L and parasitic capacitance per unit length C are 
given by [6, 7]: 

L^Lm + Lk (1) 

where Lm is the magnetic inductance per unit length and is the kinetic inductance 
per unit length of the SQUID slit given by [8]: 

K{k') 

- 2^oA^ w Aws w Aw{2w + s) 

dwk'^K^ik'Ys ""^diw + sy 2w + s d{w + s) 

and 

where //q is the permeability of free space, K{k) is the complete eUiptic integral of the 

first kind [9] with a modulus k = , k' — (1 — k'^V^^, A is the penetration depth of 

s + 2'w 

the film, sr is the dielectric constant of the STO substrate and c is the velocity of hght 
in vacuum. The SQUID slit behaves as a transmission line with distributed inductance 
L and distributed capacitance C. The impedance Zab of the slit seen from terminals A 
and B is given by [6]: 

Zab = iZo tein{niVLC) + iQLpR (5) 

where Zq — \Jl/C is the characteristic impedance of the slit transmission line, is the 
angular frequency of measurement and LpR is the junction parasitic inductance. The 
first term arises because the hairpin shaped slit can be treated as a shorted transmission 
line of length I. 
Now, using the formula [6] 



tan(f) = -^|:^^3^^^, (6) 
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eqii.(5) can be expanded as (in the lossless case with Lpji — 0): 

oo -j^ 

with Cn — CI/ 2 and L„ = 8LZ/7r^(2n — 1)^. This transformation to an equivalent cir- 
cuit has the advantage that it allows us to consider loss in the transmission line. In the 
lossy case, the rf-loss Rn is added to eqn.(7) leading to: 



where i?„ = Ln/Cn and Q is a quality factor. Here / is the SQUID slit length. This 
expression allows us to express the impedance Zab by the series of L-C-R resonant 
circuits as shown in fig. (lb). The circuit equations can be easily derived by the 
application of Kirchhoff's laws as in [1]. In fig. (lb), current entering point C should 
equal current leaving point C. We are assuming the most general case in which the 
SQUID is made up of junctions which are asymmetric [10, 12]. Here J is the circulating 
current through the SQUID inductance, Ib is the SQUID bias current, Vi and V2 are 
voltages across junctions 1 and 2 , 62 and 61 are the phases of junctions 2 and 1 and Rd 
is a damping resistance in parallel to the SQUID inductance. 
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Figure 1. (a) Geometry of the SQUID washer. Here I is the length of the sUt, w is 
the width of the electrode and s is the width of the slit, (b) Equivalent circuit of the 
washer when its parasitic capacitance distributing along the slit of the washer is taken 
into account. The circuit consists of the SQUID coupled to a series of LCR circuits. Rd 
is a damping resistance in parallel to the SQUID. Here current asymmetry is denoted 
by K, normal state resistance asymmetry by p and capacitance asymmetry by x- 

We write down the normalized circuit equations for the SQUID loop including 
random noise currents ini and i„2- Let the average junction critical current be /q, the 
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average junction normal state resistance be Rs and the average junction capacitance be 
Cj. Let the asymmetry parameter in Iq be k, that in Rs be p and that in Cj be So 
specifically, let us split [10] 2Jo = (1 + k)Io + (1 - k)Io, 2/Rs = (1 + p)/Rs + (1 - p)/^5 
and 2Cj = {1 + x)Cj + — x)Cj- We normalize currents by Iq, voltage by IqRs and 
time t by ^q/27tIqRs. The ac Josephson relation gives Vdc = dO/dr = V/IqRs , where 
is the normalized voltage and r is the normalized time. Then including random 
noise currents ini and i„2, the normalized circuit equations are: 

(1 + x)PcOi = 1{^B+ j) - (1 + pWi - (1 + At) sin(^i) - 7(^1 - 62) + ini ,(9) 

(1 -x)/3c^2 = \{^B-J) - (1 -P)^2 -(!-«;) sin(^2) +7(^1 - ^2) +^n2 ,(10) 

Here, 13c — 27rIoCjRg/^o is the SQUID McCumber parameter with $0 being the flux 
quantum and 7 = Rs/Rd- In [1] the current-voltage characteristics in the absence of 
noise was derived. It was found to be 

H = 2vdc + 7^rr^^^|(l + i^^)[dd + rf] + 2/t[df + dr] 

2|as — rrY I 

+ (1 — K^) ( [dd — rf] cos(27r0) + i[df — dr] sin(27r0) 



+ 2(5^\ds-Tf\^ ^'^ + " + + ""^ " ^"^f"" + '"^ 

- (1 - K^) - rf] cos(27r0) - i/3[sf - sr] sin(27r0)^ | . (11) 

Here d = d{Vfic)^ s = s^v^c), r = r{V(ic)-, A = A{viic) and x denotes the complex conjugate 
of X. The phase difference between the two junctions have been set to 27r0 where is 
the externally applied fiux normalized to $o- Here, 

/ = ^2vac{l + viy/^-v,, (12) 
s(u;) = (oo'^Pc + ioj) , d{uj) = A{uj) - (oo^Pc + ioJ + 2iu;-f) , (13) 
r{uj) = ipcu + xPcuJ^ , (14) 

and 

^ 1^ {-l2nI„RsVLC^%) 

Ll2nIotan{ujl27rIoRsVLC/%) 
If Lpji 7^ then A{(jj) should be replaced by A{u;)/{1 + 7rlQA{u;)Lpji/^o) in the 
calculations [1]. L and C are the SQUID inductance per unit length and SQUID 
parasitic capacitance per unit length respectively. The dielectric constant Sr enters 
through C. For more details the reader is referred to [1, 6]. 
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3. White noise 

We will now proceed to analyse the effect of noise with white power spectra, following 
closely the analysis in [2] . The calculation details are given in the appendix. The algebra 
is straightforward but very tedious. The final formula can be expressed as 

s^(fi) = ssij-lil + ai) + ss^vlil + aa) • (16) 

Here = dvdc/dih and = dvdc/dcj) and can be extracted from (11). The noise power 
per unit angular frequency are given by 

2r 

ssi,{uj)^ — , (17) 

TT 

S54>{^^) = (18) 

with r = 2nkBT/Io^Q a noise parameter, kB is the Boltzmann constant and T is the 
temperature. As in [2] we define the noise power spectra in practical units as 



= 2kBTRs 



4r2(l + «i) + ^^J(l + 27)(l + «2) 



(19) 



whose units are V'^/Hz. In the lumped case limit, the expressions in [2] are reproduced. 
The explicit formula for ai and a2 arc rather unwieldy [13]. Since we are interested in 
small asymmetry, we only quote the expressions to leading order in p, k, x with 7 = 0. 

^ cos^ 7r0 sin27r0[(ft: - p)v^cPc - pA{vdc)Vdc] (r,n.^ 

~ Srjdd Srjddss ' ^ ^ 

2sin^7r0 4sin7r0[— (k — p)f J^/?cCOS7r0 + KsssinTr^] 

" vp^dd ^ vp^ddss ' ^ 

where s,d are given in equations (13). We note that the capacitance asymmetry 
does not appear explicitly in the above formula although it is implicitly present in 
and v^. In figure 2(a), we show plots of SQUID voltage noise vs. bias voltage 
for different values of dielectric constant eji for both when the SQUID inductance 
is taken to be a lumped element and when it is taken to be a transmission line. 
The SQUID parameters arc the same as those used in [1] i.e. / = 55pm, SQUID 
inductance, Lsq{L x I) = 55pII,lQ = 4.75pA, Rs = 13.8Q, SQUID parasitic inductance 
LpB = ISpH, j3c = and damping resistance Rd = 0. Figure 2(a) shows that in 
the lumped inductance limit, the noise curves overlap with each other. Thus the noise 
as a function of bias voltage is independent of substrate dielectric constant in this 
case. However in references [1, 6, 3], it has been shown that a high substrate dielectric 
constant can cause transmission line resonances in dc SQUID characteristics. Thus, 
in case of high Tc SQUIDs, which are usually fabricated on STO substrates that are 
known to have a very high dielectric constant, it is important to model the SQUID 
inductance as a transmission line. In this case, we can see in figure 2(b), that the 
SQUID noise is definitely a function of dielectric constant and a high dielectric constant 
causes resonances to appear in the SQUID voltage noise vs. bias voltage curves at low 
voltages. This can be understood as follows. y/S^ is maximum when is maximum. It 
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Figure 2. (a) SQUID voltage noise vs. bias voltage for different values of dielectric 
constant e/j for the case when the SQUID inductance is taken to be a lumped element 
and not a transmission line. All the cases overlap which indicates that the noise value 
as a function of bias voltage is independent of en in this case, (b) Graphs for different 
values of dielectric constant cr for the case when the SQUID inductance is taken to 
be a transmission line. The dashed line is for en = 10 and the dot-dashed line is for 
cr = 100 and they overlap. The dotted line is for — 1000 and the solid line is for 
CR = 2000. 



(n — 1/4) $0 

follows from the analysis in [1] that this happens approximately at Vdc ~ 



for positive integral n. Since C is inversely related to e/j, it follows that for lower values 
for eR the first extremum in the noise occurs at higher Vdc- Resonances start appearing 
at voltage Vdc whenever the associated Josephson frequency matches the frequency of 
the lowest mode of the finite-length transmission line. 

In figure 3, we show plots of SQUID voltage noise vs. bias voltage as a function of 
asymmetries in junction parameters. Here en = 2000 is used. Again we show both the 
lumped inductance limit as well the transmission line limit. We can see that in case of a 
lumped SQUID inductance, asymmetries in junction parameters do not affect the plots 
too much and one can say that presence of asymmetries leads to a marginal increase 
in noise as a function of bias voltage. However, when we consider the transmission 
line limit, asymmetries have a significant effect on the curves especially at resonance 
positions, p asjTiimetry causes the sharpest increase. In our previous paper [1], we have 
seen that this is also the case for the dV/ dcf) vs. bias voltage curves. Thus p asymmetry 
enhances the peak in both dV/ d^) as well as voltage noise curves. 

This can be understood as follows: As we have explained in the previous section, 
the resistance asymmetry p has been chosen such that the total resistance of the 
SQUID would be a constant if the SQUID inductance behaved as a lumped element 
as opposed to a transmission line. At positions off-resonance, the conductance through 
the transmission line is aided by the finite width of the resonant peaks. Thus off 
resonance, the same physics applies and the total off-resonant conductance of the SQUID 
remains constant. From the point of view of transmission line physics, the off-resonant 
conductance is proportional to the sum of the peak widths of the arms of the SQUID. 
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Figure 3. (a) SQUID voltage noise vs. bias voltage for different asymmetry 
parameters and a fixed eu = 2000 in the case of a lumped SQUID inductance. The 
solid line is for /? = k = x = 0, the dotted line is for p = 0.2, k = x = 0, the dashed 
line is for k — 0.5, p — x — ^ and the dot-dashed hne is for k = 0.5, p = 0.2 and x = 
(b) Graphs for the transmission line case where the different types of lines stand for 
the same combination of asymmetries as in the lumped case. 



On resonance, however, the conductance is proportional to the sum of "quahty factors" , 
proportional to the inverse peak widths. This quantity increases as the asymmetry is 
increased. It needs to be emphasized that this is the asymmetry of the junctions, i.e., the 
loads of the transmission line and not asymmetries in the transmission lines themselves. 

The effect is strongest for resistance asymmetry because it enters inversely into the 
quality factor. Critical current asymmetry enters through the Josephson inductance 
(Lj (X j^). It is known that quality factor Q oc R^JC/L. Thus critical current 
asymmetry enters under the square root and thus has a much smaller influence. Also, 
both effects due to critical current asymmetry and resistance asymmetry counteract 
themselves. 

Let us now briefly discuss the noise energy defined as [2, 11] 

E = SMV2Lvl{^^)\ (22) 

Remarkably, even though this asymmetry increases the absolute noise level, the noise 
energy is generally lowered by resistance asymmetry. This is seen in figure 4 a and b. 
As introduced above, the noise energy is the appropriately normalised performance 
quantifier, relating the absolute noise to the squared transfer function. As the 
transfer function shows the same enhancement by asymmetry just discussed but enters 
quadratically into the noise energy, the resistance asymmetry reduces the noise energy 
close to resonance and thus is rather smooth across the voltages of interest. 

4. Discussion 

In this paper, we have studied the effect of asymmetry on the noise characteristics in high 
Tc dc SQUIDs which behave as transmission lines. It was shown that asymmetry can 
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Figure 4. (a) The noise energy for the SQUID plotted against various asymmetry 

parameters. The sohd line is for /? = k = x = 0, the dotted line is for p = 0.2, k = x = 
0, the dashed line is for k = 0.5, p = x = and the dot-dashed line is for k = 0.5, p = 0.2 
and X = 0- This tells us that putting in some resistance asymmetry not only minimizes 
the noise energy globally but the noise at the resonance positions is also suppressed, 
(b) Plot of noise energy vs. resistance asymmetry. This is a monotonically decreasing 
curve for p < 0.3 which tells us that noise energy can made very small by making the 
junction resistances a bit asymmetric. 

be tuned to improve the signal to noise ratio. In particular, some resistance asymmetry 
can cause a marked decrease in noise energy both globally as well as at the resonance 
positions and hence, it is not necessary to strive for extreme symmetry between the 
junctions during device design. 
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Appendix A. Calculation details 

Adding and subtracting equations (9) and (10) and denoting 9i + ^2 by 5 and 61 
we get: 

PcS" + xPcD" = iB-S' - pD - 2sin(5/2) cos(D/2) 
- 2Ksin(£>/2) cos (5/2) +813 

and 

PcD" + xPcS" = j-D'{l + 27) - pS' - 2sin(D/2) cos(5/2) 

- 2k sin(5/2) cos(L»/2) + ^(5(/) (A.2) 

Noise currents ini and z„2 have been combined to form the current noise Sib = ini + ^n2 ^^nd 
flux noise 5(p = |(ini - in2)- [2] 



- 02 by D, 
(A.l) 
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Let 5*0 and Dq satisfy equations (A.l) and (A.2) with z„i = in2 = ^ i-e- in the absence 
of noise. In the following steps, we obtain an expression for a low-frequency component of the 
voltage vs in the presence of noise where vs is defined by ^5 = We consider the case 

where the noise current Sib and the noise flux d(f) are small and express S and D as: 

S = So + Sn (A.3) 

D = Do + Dn (A.4) 

where Sn and represent variations due to the noises Sis and 5(j). 

Substituting equations (A.3) and (A.4) into equations (A.l) and (A.2), we obtain the 
following linearized equations for Sn and 

{Sn + KDn) cos(y ) cos(^) + {Dn + K^n) sin(y ) sin(^) 

= - PcS'n - xPcD'n + SlB - < - pO'^ (A.5) 



{Sn - KDn) sin( y) sin(^) - {Dn + nSn) cos(y ) cos(^) 

= PcD'n + xPcS'n + (1 + 27)D; + pS'n - |( 



(A.6) 



with 



S D 

cos( y) cos(— ^) = "^Ajn exp(zmwdci) (A.7) 

m 

So Dq 

sin( y) sin(y ) = "^3^ exp{imvdct) (A.8) 

m 

where it is assumed that S'n ^ 1 and D„ ^ 1. The trigonometric products in equations 

(A.7) and (A.8) have been expressed in a Fourier transform since the Josephson current has 

frequency components of a; = mv^c in a finite voltage state of v^c- Here, Am and Bm are 

coefficients representing the magnitude of the m*^ harmonics and ni is an integer. Since 

equations (A.5) and (A.6) arc linear, one can consider independently solutions of Sn and D„ 

for individual frequency components of of Sib and Scj). Fourier transforms of 6iB and 6(p are 

1 dS 

defined as 5Ib{'^) and 5^{uj) respectively. The voltage noise is Svg = 2~J^ ^^'^ Fourier 
transform is 6Vs{u;). Therefore, the voltage noise power spectrum is given by [2] : 

Sv{n) = (^dVs{n)dv^{n)^ (a.9) 

Next, we obtain the low-frequency component of the voltage noise, i.e. 6Vs{^) where the 
frequency a; = is considered to be much lower then the Josephson oscillation frequency 
u = Vdc- It can be shown from equations A.5 - A.8 that the low frequency voltage noise 

6Vsx{^) arises not only from low-frequency components of SIb{^) and (^$($7), but also from 
high frequency components of 6Ib{^ — mAJdc) and (^<I>(r2 — mvdc)- The voltage noise due to 
high frequency components of 51 b{^ — mvdc) and 6^{Q, — mvdc) has been expressed as the 
noise due to the Josephson mixing effect. 
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First, we obtain the value of 6Vs{i^) due to the low frequency components of 51 b{^) 
and It is difficult to solve equations A. 5 and A. 6 exactly for frequency components of 

and (^$(r2). However, since the frequency u; = i7 is much lower than the Josephson 
oscillation frequency u = ■u^o one can regard 57b (O) and (^$(0) as quasi-static changes of Ib 
and $ respectively. In this case, the value of SVs (fi) should be given by the change of the dc 
voltage Vdc due to SIb{^) and i.e. 

SVs,L{n) = ra6lB{n) + v^6^{n) (A.IO) 



where 



dvdc 



and 



dvdc 



are the dynamic resistance and the transfer function in the 



dii, d(f) 
absence of noise respectively. 

Next, we obtain the value of 5Vs{^) due to high frequency components of 5Ib and 6^^ i.e. 

5Vs^h{^) with uj = 0, — mVcic- In this regime, the R.H.S of equations (A. 5) and (A. 6) are much 

larger than the L.H.S. Thus as a first approximation to obtain the lowest order perturbation 

solution, the L.H.S are set to zero. Taking Fourier transforms of equations (A. 5) and (A. 6) 

then gives: 



6lB{oj)di 



and 



1 
P 



where, 



with 



and 



si 
di 

S2 



sl — diS2 

iujp - io^xPc 
■ + 27) - u^l3c 



(A.ll) 

(A.12) 

(A.13) 

(A.14) 
(A.15) 

(A.16) 



Equations (A.ll) and (A.12) arc the zcroth order solutions. In order to get the first order 
solutions, we plug (A.ll) and (A.12) into (A. 5) and (A. 6) L.H.S which gives expressions for 

Si and Di from (A.5), (A.6), (A.7), (A.8), (A.ll) and (A.12). Now, 6vs = ^ SVsico) = 



—ilvSniu}) 



2dt 



. Therefore, 



(A.17) 



The expression for high frequency noise with the expressions for S^ and Si substituted 
in (A.17) is quite long so we avoid writing the complete expression here. From (A.IO) and 
(A.17) we get at the measurement frequency Q, 



(A.18) 



From (A. 9), the average expectation value of (A.18) thus gives the expression for voltage 
noise power spectrum. In (A.7) and (A.8), we have taken Ai = A^i = cos(^) and 
Bi = -B_i = -isin(^) [2]. 
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